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Preface 

The  following  report  is  a  written  version  of  an  abstract  "Differential  Cartography" 
which  was  presented  on  August  21,  1991  at  the  XX  General  Assembly  of  the  Internationa! 
Association  of  Geodesy  in  Vienna,  Austria,  in  a  working  session  of  the  Special 
Commission:  Geodetic  Aspects  of  the  Law  of  the  Sea  under  the  chairmanship  of  Professor 
Petr  VaniCek  (Canada). 

The  abstract  addresses  the  various  properties  of  conformal  projections  of  a  sphere 
and  seeks  to  generalize  them  to  an  arbitrary  smooth  surface  of  revolution  by  using  the 
ba.sic  notions  of  a  new  discipline  which  we  call  differential  cartography.  In  particular,  it  is 
shown  how  the  familiar  properties  of  the  cla.ssical  Mercator  projection  may  be  carried  over 
from  a  sphere,  or  .spheroid,  to  an  arbitrary  surface  of  luvolulion. 


IV 


"DIFFERENTIAL  CARTOGRAPHY" 

Joseph  Zund 

Summary: 

The  various  properties  of  conformal  projections  of  a  sphere  and  a  surface  of 
revolution  are  derived  using  the  notions  of  differential  cartography.  It  is  shown  that 
familiar  properties  of  the  Mercator  projection  may  he  carried  over  from  a  sphere  to  an 
arbitrary  surface  of  revolution. 


1.  introducuon 

Differential  cartography  may  be  regarded  as  a  reloimulation  ol  mathematical 
cartography  from  the  viewpoint  of  differential  geometry.  Mathematical  cartography  may 
be  roughly  reduced  to  two  types  of  problems:  the 

Direct  Problem:  Construct  a  projection,  hence  a  map,  having  certain  specified  properties; 
and  the 

Indirect  Problem:  Given  certain  properties,  determine  ail  projections  having  the.se 
properties. 

In  practice,  the  Direct  Problem  has  primarily  been  the  one  which  has  attracted  the  mo.st 
attention.  Many  fascinating  projections  are  known  and  their  number  seems  limitless,  and 
restricted  only  by  the  ingenuity  of  their  discoverers.  However,  such  projections  need  not 
immediately  reveal  the  underlying  .structural  features  which  govern  their  construction. 
The  Indirect  Problem  is  concerned  with  remedying  this  situation.  Few  general  results  are 
known,  and  differential  cartography  is  intended  to  provide  a  general  approach  lor  handling 
this  problem.  Indeed,  our  primary  goal  in  the  pre.sent  paper  is  to  understand  the  Mercator 
projection,  its  principal  properties,  and  the  degree  ot  generality  intrinsic  in  its 
construction.  In  particular,  we  will  show  that  this  conlormal  projection  can  be  extended 
from  .spheres  to  arbitrary  surfaces  of  revolution  which  include  spheroids  as  a  special  ca.se. 

This  projection  was  employed  in  1569  by  Gerardus  Mercator  in  his  world  atlas. 
Nova  et  accurata  orbis  terrae  descriptio  ad  usum  navigantium  emendate  accoinmodata. 
His  world  map  measured  2x  1.32  meters,  extended  from  latitude  80°  N  to  66°  3'S, 
exhibited  parallels  and  meridians  intersecting  at  right  angles,  and  a  special  latitude  law  ol 
.scale  e.specially  devi.sed  to  yield  the  loxodrome  properly.  A  loxodrome,  or  rhumb  line,  on 
a  sphere  is  a  curve  which  intensecLs  the  meridians  at  a  con.stani  angle,  and  corresponds  to 
the  mariner's  task  of  sailing  courses  of  constant  heading.  Mercator's  problem  was  to  find  a 
projection  having  the  loxodrome  property  that  .such  courses  project  into  straight  lines  on  a 


planar  map. 

Mercator's  methodology  is  unknown,  but  it  was  undoublably  graphical  since  in  his 


lime  none  of  the  mathemuties  (e.g.  logarithms  and  the  differenlial/integral  calculus)  was 
available.  The  angle-preserving  character,  or  conformality,  of  his  projection  was  only 
gradually  recognized,  however  it  is  now  recognized  as  a  fundamental  property  of  this 
projection.  For  example,  each  loxodrome  iniensccus  every  meridian  at  a  constant  angle, 
and  under  a  conformal  projection  angles  are  pre.servcd;  it  nece.s.sarily  follows  that  the 
images  of  the  parallels  and  meridians  are  orthogonal,  and  the  image  of  a  loxodrome 
intersects  the  meridians  at  the  given  fixed  angle.  Hence,  one  may  immediately  conclude 
that  the  images  of  many  loxodromes  at  such  fixed  angles  are  parallel  straight  lines. 

Our  approach  employing  differential  cartography  is  es.sentially  a  modem  rendition 
of  classical  theory.  Lambert  (1772)  first  suggested  the  construction  of  maps  based  on  the 
use  of  the  angle-pre.serving  property.  Euler  (1777)  then  gave  the  detailed  mathematical 
theory  of  the  Mercator  projection,  and  then  proved  the  impossibility  of  obtaining  an 
isometric,  i.e.  distance-preserving,  repre.sentation  of  the  sphere  on  the  Euclidean  plane. 
Lagrange  (1779)  con.sidered  the  general  problem  of  conformal  projections  for  surfaces  of 
revolution.  The  term  ‘conformal’  was  introduced  by  Gauss  (1822,  1843)  who  studied  the 
mapping  problem  for  a  pair  of  arbitrary  .surfaces,  indicated  the  geodetic  applications,  and 
the  rich  connection  with  the  theory  of  functions  of  a  complex  variable.  The  re.sults  in 
Section  2  are  implicit  in  the  work  of  Lagrange  and  Gau.ss.  They  were  well  known  in  the 
last  century,  but  are  not  commonly  employed  in  contemporary  accounts.  Section  3 
considers  the  extension  of  the.se  results  for  surfaces  of  revolution.  1  have  not  found  them 
explicitly  stated  in  the  cla.ssical  literature,  however  they  are  in  the  spirit  of  Lagrange  and 
Gauss  and  would  have  been  no  surprise  to  them.  Our  exposition  concludes  with  a 
di.scu.ssion  in  Section  4. 

The  terminology  employed  in  t)ur  work  e.s.sentially  follows  that  employed  in 
classical  differential  geometry,  e.g.  in  the  textbook  of  Struik  (1961),  however  even  in  the 
mathematical  literature  the  u.sage  of  some  terms  is  not  uniform.  Unfortunately,  often  the 
mathematical  u.sage  of  many  icons  is  not  consi.stent  with  that  utilized  in  cartography. 
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The  following  discussion  is  given  with  some  detail  since  it  illustrates  the  basic 


approach  which  will  be  presented  in  Section  3. 

Let  S2  denote  a  unit  sphere  in  ordinary  Euclidean  3-space  E^.  We  choose  the 
Gaussian  parametrization  of  S2.  in  terms  of  the  longitude  co  and  latitude  0.  Then  the 
line  element,  or  the  first  fundamental  form,  of  S2  is 

ds^  =  co^0  dro^  +  d0^  ( 1 ) 

and  (O),  0)  determines  an  orthogonal  parameter  net  on  S2.  This  net  can  be  put  in  a 
non-trivial  isometric,  or  isothermic,  form  (see  page  171  of  Struik  (1961))  by  rewriting  (1) 
in  the  form 

ds^  =  cos^0  {dtu^  +  sec^0  d^^)  (2) 

which  amounts  to  choosing  a  unit  scale  along  the  parallels  or  a>-curves,  and  a  variable 

scale  on  the  meridians  or  0-curves  on  82- 

It  is  now  convenient  to  replace  (6),  0)  by  a  pair  of  real  isotropic  parameters 

(^,  rj).  These  are  obtained  by  employing  a  complex  factorization  of  (2): 

2  2 

ds  =  cos  0  (dtu  +  i  sec0  d0)  (dty  -  i  sec  0  d0) 
and  introducing  the  complex  variable 

C  :=  ^  +  in  •  O) 

Then 


dC  =  d^  +  idrj  =  d<0  +  i  sec0  d0 
and  so  upon  omitting  the  constants  of  integration  we  have 


(4) 

Tj  =  |sec0  d0  =  log  tan  (0/2  +  n/4) . 

2 

Thus  the  isotropic  isothermic  form  of  ds  becomes 

ds^  =  sech^Tj  +  dn^)  (5) 

since  cos0  =  sech  T], 

We  now  consider  a  conformal  projection,  or  transformation,  T  of  S2  onto  the 


Euclidean  plane  E2,  where  ’hats’  are  used  on  quantities  to  emphasize  that  they  occur  as 

A  A  A 

T-images.  Upon  choosing  Cartesian  coordinates  (x.  y)  in  E^,  the  line  element  is 
accordingly 


AO  AO  AO 

ds'^  =  dx"^  +  dy^ 


(6) 


It  is  well-known  from  complex  function  theory  that  conformal  transformations  are  defined 
by  analytic  functions,  indeed  the  conformality  of  a  transformation  may  be  taken  as 
defining  the  analyticity  of  the  function  describing  it.  Let  F  denote  an  analytic  function 
of  then  it  is  easy  to  check  that 


ds^  =  |dF|^  =  cosh^n  0) 

where  the  subscript  on  F  indicates  differentiation.  The  analyticity  of  F  is  guaranteed  by 
the  Cauchy-Riemann  equations,  which  become  F_  =  0  where  the  bar  denotes  complex 

c 

conjugation.  Since  in  differential  geometry  a  conformal  transformation  between  a  pair  of 

A 

surfaces  S  and  S  is  defined  by  their  metric  ten.sors  being  proportional,  the  same 
Gaussian  parametrization  must  he  taken  on  both  surfaces.  In  our  case  this  requires 

A  A 

identifying  (g,  rj)  iind  (x,  y)  .so  (4)  yields 


A 

X  :=  (5  =  w 

A  (X) 

y  ;=  T]  =  log  tan  {ip/2  +  n/4) . 

A  A 

The  linear  distortion  of  T  :  S2  -•  E2  is  defined  by  £j  :=  ds/ds,  and  .so  by  (7) 

£j  =  cosh  Tj  |F^|  =  sec0  |F^|  .  (9) 


The  most  familiar  kinds  of  conformal  projections  T  in  mathematical  cartography  then 
correspond  to  the  following  choices  of  the  analytic  function  F: 


a)  F  =  ^,  T  =  Tj^  {Mercator)-, 

b)  F  =  ce*^,  where  c  is  a  con.stant,  T  =  T^  {stereographic)-,  (10) 

c)  F  =  ce  ’  where  c  is  a  constant  and  k  >  0  is  a  real  number,  T  =  Tj^  {Lambert). 

Note  that  for  is  merely  the  identity  function  —  a  deceptively  simple  result  ba.sed 

on  a  rather  complicated  ^  as  in  (4)  —  and  £j  =  .sec0. 
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The  basic  results  about  Tj^  may  now  be  stated  as  follows: 

A 

(i)  Tj^  :  S2  E2  has  the  loxodrome  property  ; 

A 

(ii)  The  only  T  :  S2  -*  E2  having  the  loxodrome  property  is  T  =  Tj^. 

The  proof  of  (i)  is  immediate.  A  loxodrome  C  on  S2  intersecting  the  meridians 
at  an  angle  1?  in  the  (^,  tj)  parameters  satisfies  d^drj  =  tan  iJ ,  and  so  is  given  by  an 
isotropic  equation 

C  :  ^  =  rj  tan  d  +  c  j  ; 

which  is  eo^uivalent  to 

0)  cot  i?  =  log  tan  (<l»/2  +  7^4)  +  C2  ; 

and  whose  T-image  is 

AAA 

C  :  y  =  X  cot  t?  + 

A  A 

where  c^,  C2,  c^  are  constants.  The  equation  of  C  is  a  straight  line  in  £2-  Note  that 

A 

although  C  encircles  a  pole  of  S2  in  spiral-like  fashion  an  infinite  number  of  times,  it 
has  a  finite  length.  For  examnle  if  P  and  Q  are  the  re.spective  initial  and  final  points  on 
C,  then  the  length  of  C  is 

L(C)  =  (0Q  -  ^p)  .sec  6  . 

Hence  if  P  lies  on  the  equator  (</»p  =  0),  and  Q  is  the  (North)  pole  (0q  =  nl2)  and 

L(0  =  j  sec  j? . 

A 

Finally,  we  note  that  for  Tj.,  the  T-imagc  of  C  is  an  cqui-angular  spiral  in  E2. 

The  proof  of  (ii)  is  more  subtle  and  depends  on  the  properties  of  i.sothermic  and 
isotropic  parameter  systems.  The.se  include  the  fact  that  all  pairs  of  real  i.sothermic 
parameters  are  related  by  analytic  functions,  and  that  i.sotropic  parameters  provide  a 
particularly  convenient  parametrization  of  S2.  Indeed,  a  nece.s.sary  and  sufficient 
condition  that  a  tran.sformation  T  between  a  pair  of  .surfaces  be  conformal  is  that  their 
isotropic  curves  correspond.  This  is  the  content  of  our  equation  (8). 
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3. 


We  now  extend  the  results  of  Section  2  for  an  arbitrary  surface  of  revolution  Sj^. 
We  suppose  that  is  generated  by  the  revolution  of  the  graph  of  a  smooth  function  f 
of  a  radial  variable  about  an  axis  I  in  E^.  Choosing  I  to  be  the  z-axis,  this  radius 

2  2 

variable  becomes  J  x  +  y  and  is  denoted  by  p.  Then  the  longitude  u>  and  p  may  be 
taken  as  Gaussian  parameters  of  Sj^,  and  in  terms  of  the  parameter  net  (to,  p)  the  line 


element  of  Sj^  becomes 


ds^  =  p^do)^  +  (1  +  f'^)dp^ 


where  the  dash  denotes  differentiation  with  respect  to  p.  An  isothermic  .'orm  of  ds^  is 


given  by 


and  the  corresponding  isotropic  parameters  (4,  rj)  are 

^:=ca 


n:=|iji+r2dp. 


As  in  (3)  we  may  introduce  the  complex  variable  and  taking  g  to  be  the  inverse 
function,  p  :=  g(Tj),  of  the  second  equation  in  (13)  we  have  the  isotropic  isothermic  form 
of  ds^: 

ds'^  =  p^(d^  +  dt]") 

=  g^(d5^  4  dt)2)  (14) 

A 

Let  T  :  Sj^  -4  E2  denote  a  conformal  transformation,  so  as  in  Section  2  we  have 

^^9  I  99 

ds  =-4|F.1V  (15) 


where 


A 

x:=^ 


y  :=  n- 


fi 


The  linear  distortion  of  T  is  then 


(17) 


Note  that  in  (17)  the  function  g  characterize  the  particular  form  of  Sj^,  and  F  is  an 
arbitrary  analytic  function  of 

The  analysis  in  Section  2,  including  the  specializations  of  F,  indicated  in  (10) 
carry  over  to  define  and  Tj^  respectively.  Hence,  we  immediately  have  the 

basic  results: 


(i)  Tj^  :  -t  E2  has  the  loxodrome*  property, 

A 

(ii)  The  only  T  :  Sj^  -*  E2  having  the  loxodrome  property  is  T  =  Tj^. 

The  proofs  of  (i)  and  (ii)  are  completely  analogous  to  the  reasoning  given  in  Section  2 
when  (4)  is  replaced  by  (13).  It  would  be  truly  remarkable  that  the  above  results  were  not 
established  by  the  mathematical  cartographers  of  the  nineteenth  century.  However,  we 
have  been  unable  to  locate  specific  references  in  the  literature  which  e.stablish  this.  In  any 
case,  our  approach  which  is  based  on  differential  cartography  and  properties  of  analytic 
functions  may  well  be  new. 


4.  Discussion 

The  results  in  Section  3  may  be  surprising,  however,  and  if  this  is  so,  it  is  bccau.se 
the  familiar  results  on  the  Mercator  projection  in  Section  2  are  not  usually  di.scussed  as  a 
Indirect  Problem.  The  viewpoint  of  differential  cartography,  indeed  la  raison  d'etre  for  it, 
is  to  provide  a  geometric  approach  to  mathematical  cartography.  In  this  rc.spect,  the  key 
ingredients  leading  to  the  results  (i)  and  (ii)  arc  the  emphasis  on  the  i.sothermal  form  and 
isotropic  parameters  for  the  surface  line  element.  These  notions  are  purely  of  a  geometric 
character  and  their  full  power,  e.g.  their  connection  with  complex  function  theory,  would 
seem  to  be  evident  only  within  the  framework  offered  by  differential  cartography. 


•The  notion  of  a  loxodrome  can  be  generalized  for  an  Sj^  to  be  a  curve  which 
intersects  a  pencil  of  planes  at  a  constant  angle.  The  axis  I  of  Sj^  is  then  taken 
as  the  axis  of  the  pencil  which  ct)n.si.sts  of  meridional  planes. 
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Three  conclusions  are  immediately  suggested  by  the  analysis  given  in  Sections  2 
and  3.  First,  by  viewing  the  conformal  projection  as  being  defined  by  an  analytic  function 
F  of  a  complex  variable,  it  is  clear  that  the  three  cases  indicated  in  (10)  represent  only  the 
simplest  choices  of  F.  Indeed,  on  function-theoretic  grounds  it  could  be  argued  that  these 
choices  are  almost  the  simplest  possible  analytic  functions.  Hence,  the  Mercator, 
stereographic,  and  Lambert  projections  by  no  means  exhaust  the  possibilities.  Thus,  the 
use  of  conformal  projections  in  mathematical  cartography  is  hardly  a  fully  explored 
technique,  until  other  choices  of  F  are  investigated. 

The  second  conclusion  is  that  the  results  for  Sj^  offer  the  possibility  of  replacing 
the  spheroid  by  a  pear-shaped  globe  having  a  hump  at  the  North  pole  and  a  depression  at 
the  South  pole  provided  that  these  irregularities  can  be  regarded  as  being  axially 
symmetric.  In  any  ca.se,  one  could  devi.se  a  conformal  projection  which  could  be  more 
accurate  than  the  usual  spheroidal  model.  Indeed,  one  could  employ  Tj^  in  such  a  case 
by  taking  the  base  parallel  not  as  the  equator  (0  =  0),  but  as  a  higher  latitude  (f>  =  tp^\ 
Recently  VaniCek  and  Sjdberg  (1991)  have  considered  an  analogous  generalization  of 
zonal  harmonics  in  seeking  a  new  description  of  the  gravity  field  of  the  Earth. 

Finally,  (ii)  of  Sections  2  and  3  .show  that  for  conformal  transfonnations  Tj^  is 
the  only  projection  having  the  loxodrome  property.  However,  this  property  requires  that 

A 

angle  t?  be  preserved  under  a  mapping  of  S2,  or  Sj^,  into  E2,  hence  it  is  meaningful 
only  for  a  conformal  transformation.  Thus,  we  may  conclude  that  the  loxodrome  property 
can  hold  only  for  a  conformal  transformation  T,  and  as  already  noted  by  our  item  (ii)  we 
must  have  T  =  Tj^.  In  i)ther  words,  the  loxodrome  property  uniquely  characterizes 
Mercator  transformations. 

I  am  indebted  to  Profe.s.sor  Petr  Vanicek  for  his  suggestion  to  investigate  this 
question,  and  to  Bernard  Chovitz  for  his  caieful  .scrutiny  of  a  preliminary  version  of  the 
manuscript  and  several  u.seful  di.scu.ssions. 
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